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Abstract 

We resum Sudakov threshold enhancements in heavy quark hadroproduction for single- 
heavy quark inclusive and pair-inclusive kinematics. We expand these resummed results 
and derive analytical finite-order cross sections through next-to- next-to-leading order. 
This involves the construction of next-to-leading order matching conditions in color space. 
For the scale dependent terms we derive exact results using renormalization group meth- 
ods. We study the effects of scale variations, scheme and kinematics choice on the partonic 
and hadronic cross sections, and provide estimates for top and bottom quark production 
cross sections. 



1 Introduction 



Long- and short-distance dynamics in inclusive hadronic hard-scattering cross sections are fac- 
torized in Quantum Chromodynamics (QCD) into universal, non-perturbative parton distri- 
bution functions and fragmentation functions, and perturbatively calculable hard scattering 
functions. Remnants of long-distance dynamics in a hard scattering function can, however, 
become large in regions of phase space near partonic threshold and dominate higher order 
corrections. Such Sudakov corrections assume the form of distributions that are singular at 
partonic threshold. Threshold resummation organizes these double-logarithmic corrections to 
all orders, thereby extending the predictive power of QCD to these phase space regions. 

Early on [I], 0], the organization of such corrections to arbitrary logarithmic accuracy was 
achieved for the Drell-Yan cross section. An equivalent level of understanding for general 
QCD processes with more complex color structures at the Born level has been achieved more 
recently || £|, ^ ||, [F], || . The resummation of Sudakov corrections in such processes to next-to- 
leading logarithmic (NLL) accuracy requires understanding how these structures mix under soft 
gluon radiation. Many NLL-resummed cross sections have been calculated: heavy quark hadro- 
0- II- fl Mi OH 0' HH an d electroproduction |l3|, [yj, dijet production [|, f?j, single-jet production 



15f , Higgs production [|nj , prompt photon production [TJ, |17|, [L8|, [19], |20| and hadroproduction 
of electroweak bosons ||21[| . For a recent review see Ref. |p2| . The formalism of Refs. || |], ||, |^, [7|] 
and the oneQ of Ref. [Q both allow arbitrary logarithmic accuracy. Processes involving Born- 
level two-particle scattering may be described in either single-particle inclusive or pair-inclusive 
kinematics. NLL resummation was initially performed in pair-inclusive kinematics 0, [S], |6], [7| 



and later extended to single-particle kinematics [12 , 17 1. 

Resummed cross sections constitute an approximate sum of the complete perturbative ex- 
pansion if, at each order, the Sudakov corrections dominate. The numerical evaluation of 
resummed cross sections requires a prescription to handle infrared renormalon singularities. 
Resummed cross sections may also be expanded to provide estimates of finite higher order 
corrections which do not suffer from renormalon problems. In this paper we shall employ 
the resummed cross sections in the latter fashion: as generating functionals of approximate 
perturbation theory. 

Resummed results for heavy quark production at leading logarithmic accuracy have been 
presented some time ago p4|, |25], |26|, [27|, |2^, p9|| . In our paper we expand the NLL resummed 
cross sections presented in Refs. [B, [], [E| and derive complete analytic expressions through 
next-to-next-to-leading order (NNLO) for double-differential heavy quark cross sections in two 
different kinematics: heavy quark pair-inclusive and single-heavy quark inclusive. To achieve 
next-to-next-to-leading logarithmic (NNLL) accuracy we include color- coherence effects and 
contributions due to soft radiation from one-loop virtual graphs via matching conditions. 

Our paper is organized as follows. In section 2 we discuss both types of kinematics and 
their singular functions. Section 3 describes the construction of the resummed cross sections. 
In section 4 we expand the resummed cross section to NLO and NNLO and present analytical 
NNLL double-differential cross sections at each order, with both types of kinematics, in the 
gluon-gluon (gg) and quark- ant iquark (qq) channels. We numerically study the inclusive par- 
tonic and hadronic cross sections in sections 5 and 6 respectively. Our conclusions are presented 
in section 7. Appendix A contains the NLO matching terms for pair-inclusive kinematics. In 
appendix B we collect all our explicit expressions for the NLO and NNLO differential cross 



1 A recent study, Ref. f23fl, compares these formalisms for prompt photon production. 
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sections. 

A companion study |3(|, also addressing third and fourth order corrections, contributions 
from subleading logarithms, and some differential distributions, already featuring some of the 
second order results for inclusive cross sections derived in this paper, was recently presented by 
one of us. 



2 Kinematics and cross sections 

The kinematics of inclusive heavy quark hadroproduction depend on which final state momenta 
are reconstructed. In threshold resummation this kinematics choice manifests itself at next-to- 
leading logarithmic level |TI| . We discuss two types of near-elastic kinematics in heavy quark 



hadroproduction, one-particle inclusive (1PI) and pair-invariant mass (PIM) kinematics. 

2.1 One-particle inclusive (1PI) kinematics 

Heavy quark hadroproduction in 1PI kinematics is defined by 

hx(Px) + h 2 (P 2 ) — > Q(pi)+ X[Q]( Px ), (1) 

where hi and h 2 are hadrons, X[Q] denotes any allowed hadronic final state containing at 
least the heavy antiquark, and Q(pi) is the identified heavy quark with mass m. The hadronic 
invariants in this reaction are 

S = (P 1 + P 2 ) 2 , T 1 = (P 2 -j9 1 ) 2 -m 2 , Ul = (P 1 -p 1 f-m 2 , (2) 

and 

S 4 = S + Ti + U x , (3) 

where S4 is a measure of the inelasticity of the hadronic reaction ([]]). Near threshold, reaction 
([!]) is dominated by the partonic subprocesses 

q{kx)+q{k 2 ) — > Q(pi)+X'[Q](p' 2 ), (4) 
g(kx)+g(k 2 ) — Q(pi)+X'[Q](p 2 ). (5) 

If X'[Q](jo 2 ) = Q,(p 2 ), the reaction is at partonic threshold and the heavy antiquark has mo- 
mentum p 2 . Note that threshold production does not mean that the heavy quarks are produced 
at rest. The qg and qg channels contribute at one order higher in a s than the reactions (|J) and 
The partonic invariants corresponding to (@) are 

s = (k 1 + k 2 ) 2 , tx = {k 2 - pi) 2 - m 2 , u x = (h - pi) 2 - m 2 . (6) 

The invariant S4 = {p 2 ) 2 — m 2 which measures the inelasticity of the partonic reactions (|J) and 
d^) is related to the other partonic invariants by 

s 4 = s + tx + Ux- (7) 



2 For example, for top quark pair production at the Tevatron pp collider these subprocesses contribute only 
1% of the total cross section pi], 51, |33|. 
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The inclusive partonic cross section may be calculated from 



cr t ,(s,m 2 ) = I di-h) I ds 4 ggg^lil) ( 8 ) 



s(l+/3)/2 
s(l-/3)/2 

where (3 = \J\ — 4m 2 / s and 

2 

s - ax = s + t 1 + ^. (9) 



The recoil momentum p' 2 may be split into the momentum at threshold, p 2 , and the momentum 
of any additional radiation above threshold, k, i.e. p' 2 = p 2 + k. Then, when k 2 is small, we 
can define a dimensionless weight wipi |3[ that measures the distance from threshold in 1PI 
kinematics which can in turn be expressed in terms of a vector £f PI : 

s 4 2p 2 -k 2CiPi-fc , 1fV . 
iwipi = —= ~ ^— = . (10) 



Ref. 33 contains an exact NLO treatment of this kinematics at the parton and hadron levels. 



2.2 Pair- invariant mass (PIM) kinematics 

Heavy quark hadropro duct ion in pair-invariant mass kinematics is defined by 

/ia(Px) + h 2 (P 2 ) — QQ(p') + X{p x ) . (11) 

At the parton level, the important reactions are 

q(h)+q(k 2 ) — > QQ(p') + X'(fc) , (12) 
+ g(k 2 ) — > QQ(p') + X'(fc) , (13) 

where p' 2 = M 2 is the pair-mass squared. If X'(k) = 0, the reaction is at partonic threshold 
with M 2 = s. Then 

M 2 

h = — — (1 - Pm cos9) , 
M 2 

ux = — — (l + /? M cos0) (14) 



where (3m — y 1 — Am? /M 2 and 6 1 is the scattering angle in the parton center-of-mass frame. 
In PIM kinematics, the inclusive partonic cross section may be calculated from 



S[1+ r )/2 r d 2 rr fs M 2 t ) 

a^m 2 ) = / d(- tl ) / dM 2 da t S ;Z ■ (15) 



s(l-/3)/2 

with 



The weight u>pim measures the inelasticity in PIM kinematics, 



s-M 2 2p'-k 2( P1M -k n ^ 

wpim = 1- z = ~ = -= — , (17) 

s s V s 

where z = M 2 /s. Equation ([[7]) defines the vector Cpim i n terms of the heavy quark pair 
momentum at threshold (indicated by the bar). An exact NLO treatment of this kinematics at 



the parton and hadron levels may be found in Ref. |p3 | 



2.3 Inclusive cross section 

The inclusive hadronic cross section is obtained by convoluting the inclusive partonic cross 
sections, Eqs. (||) and (|15|), with a parton flux factor 

i i 

.2\ _ _ f J„ f j„ IV _\ i . /„ ..2\ / / . 2^ 



$ij(r,fi ) = t J dxi J dx 2 8(x 1 x 2 -r) c^i/h^X!,^ )(j)j/ h2 (x 2 ,fi ) , (18) 



where <fri/h{x,fi 2 ) is the density of partons of flavor % in hadron h carrying a fraction x of the 
initial hadron h momentum, at factorization scale /i. Then 

a hlh2 (S,m 2 ) = J — ^ij( T ^ 2 ) cr lj (TS,m 2 1 y 2 ) (19) 



E / rflogio^^— ln(10) ^(r/,/! 2 ) ay(r/,m 2 ,/i 2 ) 

- „ J —oo 1 + Tj 



where 

s tS , . 

The sum is over all massless parton flavors. The second equality in Eq. ([B5]) facilitates inter- 
pretation of the figures in section 5. 

2.4 Singular functions 

The double-differential hadronic cross section for reactions (|T|) and (|TT|), d 2 a hlh2 /dV dW , enjoys 



the factorization j35| 



d 2 <Jh 1 h 2 {S, V, W) ^ f dx\ r dx 2 



5 = ^ j — j—<t>i /h M^)<t>HH^2^) 

x Uij{w K , s,t 1 ,ui,m 2 ,fj, 2 ,a s (fj,)) + O (J^j > (21) 

where is the partonic cross section, or hard scattering function, whose dependence on 
kinematics is indicated by the weight wk with K = (1PI,PIM). The variables V and W 
represent either 1PI or PIM kinematic variables such as the transverse momentum and rapidity 
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of either a single heavy quark or the heavy quark pair, respectively. The parton momentum 
fractions x\ and x 2 have lower limits x~{ and x^ which depend^ on V and W. 

We shall resum the higher order logarithmic contributions to Uij that are singular at thresh- 
old. The arguments of these logarithms are the weights Thus, the 1PI singular functions 
are plus distributions in S4, 



ln z (s 4 /m 2 
s 4 



hm { — — J -6(s 4 - A) 



1 



l + l 



In 



1+1 



(4) 



(22) 



while the PIM singular functions are plus distributions in 1 — z 

V(l-z) 



~\n\l-z) 


= lim < 


1 - z 





1-2 



-0(1-2-5) + 



1 

ITT 



ln ,+1 (5) <5(1 



(23) 



These singular functions yield finite large results when convoluted with smooth but rapidly 
changing functions such as parton densities. Note that we have normalized the 1PI functions 
to have mass dimension —2. Because our paper deals mostly with finite order cross sections, 
we denote corrections as leading logarithmic (LL) if I = 2i + 1 at order 0(a l s +3 ), i — 0, 1, . . ., 
as next-to-leading logarithmic (NLL) if I = 2i, etc. 

It is often convenient to work in moment space, defined by the Laplace transform with 
respect to wk 



f(N)= I dw K e- NwK f{w K ) 



(24) 



The upper limit of this integral is not very important, and may be set to 1, where Iuwk = 0. 
Under Laplace transformations, the plus distributions in Eqs. (p2|) and (^) become linear 
combinations of ln fc (iV) with k < I + 1 and N = N exp(7#) where 7^ is the Euler constant^. 
The precise correspondence to second order can be found in Ref. @, and through fourth order 
in Ref. []30]]. We shall almost always talk about logarithmic accuracy in momentum space. 
Occasionally in the following we denote in moment space leading logarithmic corrections at 



order Oia 



i+2\ 



1,2, .. . to correspond to I = 2i, NLL ones to I = 2i — 1 etc. Although 



the Laplace transformation of the Z'th plus distribution generates lower powers of In iV besides 
\n l+1 N, there should be no confusion in practice. 

We work in axial gauge, n • A = 0, with the gauge vector chosen as n M = (£. The implicit 



renormalization scheme is that of Ref. |$6|] and m is a pole mass. The renormalization scale hr 
is assumed to be equal to the factorization scale /x, except where explicitly indicated otherwise. 
The scale dependence of the coupling constant is controlled by the QCD /3-function 



dfl 7T 



(5(a s (jj)) 



-2bo 



26, 



71 



71 



(25) 



where b 2 = (11C A ~ 2n/)/12, and 63 = (34C^ - 2n f (3C F + 5C A ))/48. 

3 When V, W = Ti, Ui in 1PI kinematics one has x± = -Ui/(S + T t ) and = -x i T 1 /(xiS + Ui) while for 
V, W = M 2 ,Y in PIM kinematics one has x 1 = M exp(Y) / V~S and x 2 — M exp(— Y)/\/~S where Y is the pair 
rapidity. 

4 A11 the TV's in the remainder of this work are actually N, unless specified otherwise. 
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3 Resummed cross sections 



Here we describe the threshold resummation of the heavy quark hadroproduction cross section 
in PIM P, m and IPI [|T!J kinematics. Both resummations can be presented simultaneously 
using the methods and results of Refs. || [|, ||, ||, 0, |I2|| . 



3.1 Refactorization 

The resummation of the singular functions in Eqs. (p2|) and (|23|) in the perturbative expansion 
of Uij rests upon the refactorization of cjjj into separate functions of the jet-like, soft, and 
off-shell quanta that contribute to its quantum corrections. This refactorization, valid in the 
threshold region of phase space, is pictured in Fig. [l|. Each of the functions ip, h, and S organizes 




Figure 1: Refactorized form of heavy quark partonic cross section near threshold. 



large corrections corresponding to a particular region of phase space. The meaning of the /, J 
indices will be given shortly and the coefficients c^-j are given below in Eq. (|43|) . Factorizations 
of this type have been discussed earlier for deep-inelastic scattering and Drell-Yan production 
|], heavy quark hadro- and electroproduction j|, H dijet @ and Higgs |16| production, 



and single-particle/jet inclusive cross sections ^TJ. They may be generalized to include 
recoil effects [ j37j] . 

Figure [l] indicates that each factorized function depends on its own weight function ||. In 
essence, the choice of working in IPI or PIM kinematics depends on how the total weight wk 



is constructed from the individual contributions of each weight function [y, |12|1 : 



PIM : Wpim — Wi + Wj + w s (26) 

IPI : w 1PI = Wl ( ^) + Wj + w s . (27) 



In the case of IPI kinematics it is convenient to define 



IPI: N U = N(-^), N t = N f — ^ ) • (28) 



9 / * 1 \ 9 
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In both kinematics, the moments of the ij partonic cross section, Eqs. (Q), @, (|T2|), and (p~3|) 
can be written in refactorized form, up to 0{1/N) corrections, as P, j5[, |6|, |l"2" 



uj ij (N,s } t u ui } m 2 ,ij, 2 } a s (iJ,)) = h*j j(( K , m 2 , /j 2 ) ^'^^~>C^j hjiibc, m 2 , /j 2 ) 

j>i/i(N u , kj. ■ Ck/v) ^j/j( N t, h ■ Ck/p) ^ 

The indices^ I and J take values in the color-tensor space spanned by the invariant tensors 
that combine the SU(3) representations of the external partons at threshold into a singlet. The 
vector defines the kinematics, see section 2. The "incoming-jet functions" ipi/i describe the 
dynamics of partons moving collinearly to the incoming parton i. The distributions (piu are 
defined at fixed light-cone momentum while the functions ip^i are defined at fixed ki ■ (k- The 
(real) ipiu functions include all leading and some next-to-leading singular functions and are 
diagonal in color-tensor space. The function Sijji, a Hermitian matrix in color-tensor space, 
summarizes the dynamics of soft gluons that are not collinear to the incoming partons and 
contains the remaining next-to-leading contributions. Note that (next-to-leading) singularities 
associated with soft radiation from the outgoing heavy quarks^ are included in Sijjj @, ||. 
Finally, the Hermitian matrix H^u = hijih^j incorporates the effects of far off-shell partons 
and contains no singular functions. 

The jet and soft functions in Eq. ( p9|) can each be represented as operator matrix elements 
jl, ||, H B HI]- The refactorization of the cross section in Eq. fl29|) can in fact be seen as a 
separation of near-threshold degrees of freedom into distinct effective field theories J5L [|, 0]. In 
each effective theory a resummation of the singular functions may be performed via appropriate 
evolution equations. The resummed cross section is composed of contributions from all these 
effective theories which must be matched together properly, at a specified scale and to a certain 
order in perturbation theory. In this paper we match to NLO, leading to NNLL accuracy in 
our finite order expansions. Our matching scale is the heavy quark mass m. The matching 
procedure is described in section 4. 

In the following we describe the resummation of the jet and soft functions. Most of these 
steps have been discussed previously in the literature 0, ||, ||, [12| so that our description is 
brief. 



3.2 Resummed Jet and Soft Functions 

The (exponentiated) N dependence of the ratio ifii/i/fa/i in Eq. ( p9| ) follows from its factorization 
properties [|. The function ipi/i obeys, beside a renormalization group equation, an evolution 
equation governing the energy dependence, expressed as gauge-dependence |], || [38|, |39|. Solving 



these equations leads to the resummed form of this ratio. The function ip has been defined in 
PIM kinematics for i = q [J] and g || with the 1PI equivalents given in Refs. |13|, |12 |. The 



resummation of Sudakov logarithms in 1PI kinematics was verified [|13j to trace those in PIM 
kinematics 0. Expressed in terms of Wi, the ipiu are the same in PIM and 1PI kinematics. 

5 Unless stated otherwise, repeated indices are summed over. 

6 Were we to treat the heavy quarks as massless, each heavy quark would be assigned its own jet function 
|], 0, 0. However, because the heavy quark mass prevents collinear singularities, all singular functions 
arising from the final state heavy quarks are due to soft gluons. 
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The resummed ratio if>i/i/(f>i/i is 



In 



(N, 2kt ■ Cjc///) 



H=2ki-C, K 



Et (N, 2h ■ (k) 



with the MS exponent 

Ei (N, 2h ■ Ck) 



dw 



'l-e 



-Nw\ 



w 



f d\ 

J —A t [a s (V\2ki-( K ) 



(30) 



(31) 



+ -v i [a s {w2k i -( K )\ . 



The functions A\ and v l differ by color factors for i = q and g so that 



a. 



71 
71 



a. 



71 



—C F +[-) [~C F K)+..., A g (a s ) = —C A + (-) [-C A K)+... 



(2C F ) + ... 



71 



71 



a, 



71 



{2C A ) + ... 



(32) 



with K = C A (67/18 - tt 2 /6) - 5n f /9 

To incorporate the effects of scale changes in finite-order a s expansions, we note that the 
ratio ipi/i/ 4>i/i transforms under renormalization scale, //r, and factorization scale, //, as 



d . fa/i (N, 2k { ■ (k/hr) 

jiR- — In — = 

dfi R (t>i/i(N,n) 

d ipi/i (N,2ki ■ (k/i^r) 



\i— In 

djj, 



<f>i/i(N,fi) 



7* = 2-fi (a s (fx R )) 
-2%/i (N,a s (fx)) , 



(33) 
(34) 



where 7$ and 27^ are, respectively, the anomalous dimension of the quantum field i and of 
the operator whose matrix element represents HI], [12] the MS-density The anomalous 

dimension for the renormalization scale dependence of [[J in Eq. ([£J) does not depend on 
the moment N because the ultraviolet (UV) divergences of ij)iu are only due to wavefunction 
renormalization of parton i. Hence 7^ = 2 / ji where 7, is calculated in the axial gauge. The 
anomalous dimension 7^ controlling the factorization scale dependence does depend on N. 
The axial gauge anomalous dimensions are (neglecting 0(1/ N) terms) 



lq M 
Iq/q {N, a s ) 

lg («») 

Ig/g (N, a s ) 



3 n a s 
-U F h 

4 71 



-\C F \nN - - C F 



71 

b 2 h • • ■ 



71 
71 



C A lnN-b 2 



71 



1 



-C A K In N) + 



(35) 

(36) 
(37) 

(38) 



where 62 is given below Eq. ( |25| ) and K is given below Eq. (|32|). 
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The composite operator |7j that defines Sijjj contains UV divergences beyond the self 
energies of its external legs. An additional renormalization removes these extra divergences. 
However, the product Hiju§ijjj has no UV divergencies beyond those taken into account 
by the renormalization of the elementary fields and couplings of the theory, because the 
renormalizat ion-group invariant cross section and this product only differ by the functions 
which are renormalized as if they were (the "square" of) elementary fields (|33p. Hence the extra 
UV divergencies of Sijji are balanced by similar ones in H^u. Therefore, the unrenormalized 
expressions Sf^jj and H^jfJ renormalize multiplicatively ||, |^. |(| 

H i*ij = Zr l Zj l [z ijS ) IK H i jKL{zl jS ) Lj , (39) 
S$?i = (4s) JK S ijKL (z ijS ) Li . (40) 

The factors Zj denote the renormalization constants of the external fields i = q,q,g. For a 
given parton pair ij, the (Zijs)u constitute a matrix, in color-tensor space, of renormalization 
constants for the overall renormalization of the soft function. Note that (Zijs)u also includes 
the wave function renormalization of the external heavy quark legs. 

Equation P0| ) leads directly to a renormalization group equation for the matrix S^jj ||, |]. |J 



d 

l dfi R 



^R-TTT J I - ~F%JK $ij KI - Sij JK T^ KI . (41) 



The soft anomalous dimension matrix T l g U is obtained from the renormalization constants 



(Zs)u, computed in d = 4 — e dimensions, by |43 

9s & 

T hj = - y ■7^-Res e ^o{Z ijS ) IJ {a s , e) . (42) 

with g 2 s = 4na s . The solution of Eq. (|4lD is in general expressed in terms of path-ordered 
exponentials 0, |£J , see below. 

An explicit expression for 7J requires a choice of basis tensors in color-tensor space. We 
choose an s-channel singlet-octet basis 

\pqql)mnkl = {$mn$kl j {Tp) n m (Tp)k,{) , 
^ggl)abkl (fiabfikl , d a bc(Tp)kl ,ifabc(Tp)kl) , (43) 

where T F are the SU(3) generators in the fundamental representation, where the indices 
m, n, k, I take values. The indices a, b, c are adjoint indices while d a b c and f a b c are the to- 
tally symmetric and antisymmetric SU(3) invariant tensors, respectively. In this basis, the 
one-loop soft anomalous dimension matrix for the qq process is [f|, ^. ^] 



1 s 



-pqq -pqq 

1 11 1 12 

-pqq -pqq 

1 21 1 22 



(44) 



with matrix elements 



Til = s -C F [Lp + ln(2^^)+m] 

71 

rgf = 2 -^ln(^ 

71 \ti 
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1 22 



a s Cf 
~~C~a 

a -\c F 

7T I 



4 In ( y) ~ H^V^) ~ Lp - m 



+ 



c. 



-3 In 



-In 



' m 2 s 

JlUi 



+ L g + ni 



The function Lg is denned as 



The variables Ui, i — 1, 2 are 



1 - 2m 2 /s (. 1-P 







V; 



2(k t • C 



1 + P 



+ 7TZ 



(45) 



(46) 



(47) 



The one-loop soft anomalous dimension matrix for the gg process [[5], ^2) can be written as 

rff o 

^99 



-p99 
1 S 



2 1 31 

V c p39 
X 1 31 



~p99 N^-i ^gg r fi^ 
1 31 ~WT 31 1 22 



(48) 



where the three independent matrix elements are 



rff = 


a s 


7T 


rff = 


2a 




r99 _ 

1 22 — 


a s 




7i 



-C F (Lp + 1) + C A (-- In (4uiu 2 ) + 1 - tti 



(49) 



In 



-C F {Lp + 1) + ^ 



m 2 s 



+ Lp — ni 



-In (4z/!Z/ 2 ) + 2 + ln 

Note that we did not absorb the function v % of Eq. (|3lf) into Tg. 

Combining the solution of Eq. (PO) with the resummed jet functions and the hard function 
and using matrix notation for and SV,-, the full resummed partonic cross section is 

S (N, s,t 1 ,ui,m 2 ,fi 2 ,a 8 (fi)) = Trj m 2 , m 



-res / 



(50) 



xP exp 



/i' 



5^(1, Or) P exp 



exp (Ei(N U} fj,, fj, R fj exp (Ej(N u fj,, /i R )J exp I 2 J — (a s (//)) + 7j (a s (// 



i/iV (7,,' . . 

/i' 

d// 
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where the trace is in color-tensor space, and P refers to path-ordering in jji! such that (a s (m)) 
is ordered to the far right while T l g(a s (m/N)) is ordered to the far left. The operation P orders 
in the opposite way. The exponential exp(Ei) is given by 



exj>(Ei(N u , fi, fi R )) = exp 



Ei(N u ,2ki-C K 



(51) 



x 
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The exponential exp(Ej) is identical except for obvious relabelling. In Eq. ( pO]) we used the 
renormalization group behaviour, derived from Eqs. (|39|) and fl4"0|), of the product HijjjSijji 



d 



dfx 



In 



#tj / J (Vr) Sij J I (N,/J,r) 



-27< (a 5 (^ B )) - 27i (a s (///i)) 



(52) 



where we have only indicated the \ir and iV dependence. Notice that all dependence cancels 
in Eq. @). 

So far, the results have all been given in the MS-scheme. The qq results can also be presented 
in the DIS scheme, in which the function E q of Eq. (|3lD is 



E q (N, 2ki-( K ) 



E q (N, 2h ■ (k) 



DIS 



MS 



+ dw ■ 



1-e 



d\ 

T 9 



a s (vAs) 



+ <7, 



DIS 



a s (-\Ms) 



(53) 



with 



^ DIS 



(««) = -~ A C F— + 

4 7T 



(54) 



In some cases, e.g. for the unexpanded resummed cross section, it is preferable to diagonalize 
the matrix Y l g u []7], [TT], ^ [0] by a change of basis in color-tensor space, 

T^s = R- 1 T s R, (55) 

written in matrix notation with T^g diagonal. The channel labels have been suppressed. The 
columns of the matrix R are the eigenvectors of IV This leads to ordinary exponentials in 
the solution of Eq. ([£!]), rather than path-ordered ones. The elements on the diagonal are the 
(complex) eigenvalues A/. In the new basis, indicated by the subscript R, the matrices S and 
H are 



S f 



R SR , Hp = R ^ H ( R 1 



again in matrix notation. Eq. ( fH"|) then becomes 

d ~ 



Sr,ji — — v^i + Kj) Sr,ji 
djj, 



(56) 



(57) 



The matrix R qq is given in Refs. pTL E2[ while R gg appears in Ref. [22]. The solution of Eq. fl5 



is 



m 



Sr,jA-^-,CkJ = S RjJI (l,( K ) exp 



m/N 



^{A / (a s ( / i')) + A}(a s (/i'))} 



so that (reinstating the channel labels) 



m ,m 



(58) 



(59) 



i, j 



x 



SR,ijji(l,(K) exp{% J7 } exp 2 J -7- {7* (««(a* 7 )) +7i ("-(^O)} 
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The exponential in Eq. fl59|) carrying color-tensor indices is 

exp^j/) = exp(Ei(N u , fx, jj R )) exp(Ej(N t , jj, n R )) (60) 

m/N 



x exp 



/ -^{A / («,(/x')) + A}(a s (/i / 



with exp(_Ej) and exp(Ej) given in Eq. (|51~D. We derived the bulk of our subsequent results in 
both bases. 

4 NNLL-NNLO expansions for partonic cross sections 

In this section we derive analytical expressions for partonic double-differential heavy quark 
cross sections up to NNLO by expanding the resummed cross section of the previous section. 
We concentrate here on the derivation; our explicit results are collected in the appendices. We 
obtain expressions for both 1PI and PIM kinematics and in both MS and DIS factorization 
schemes. Our aim is NNLL accuracy, as defined in section |2.4j, including the scale dependent 



sector where we distinguish between the renormalization scale [Ir and mass factorization scale 
fi. This means that for coefficients of ln J (/i 2 /m 2 ) or In* (fi R /m 2 ) we include the most singular 
plus distribution and the two next-most-singular ones. 

To achieve NNLL accuracy we must derive NLO matching conditions for certain functions 
in the resummed cross section, Eq. (|5(J). To this end we expand the factors in this equation and 
identify order by order the perturbative coefficients of the functions in Eq. ([)(]). The coefficients 
can in part be explicitly computed and in part inferred by matching to exact NLO results. 

4.1 Expansions 

Let us expand each factor in Eq. ( pU[ ) in powers of a s . We set fiR = m here because Eq. ( |5UD is 
independent of {ir, leading to a cross section expansion in a s {m) which can easily be changed 
back to a s (fiR). Neglecting 1/N terms, the two- loop expansion of exp(Ei) may be written as 

exp^iV^m)) ~ 1 + ^ (j:ci: W \n k (N u )] + f^V ( £ C£ (2) \n k (N u )) + . . . (61) 



n \k=o I V7r ' \fe=0 



The coefficients C 1 ^ can be computed from the results in section [| To NNLL accuracy they 
are 

r i,{l) _ Al) r i,(l) _ 1 „ , A (l)j r i,(l) _ t r i,(l) 



CP = 4« + c 2 (4 1, ) 2 ^-5('= i ) 2 + ^ + iK , ) 2 i^ 

= (cAA^-ld^ + A^i. + i-h^-d^Aii 



2"' 1 ' 2" 



Co^ = (h 2 d i+ ld*)ll, (62) 
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where we have used the notation 



u i,(D_ 2^.(D l n 



m 



d„ 



d g = b 2 , 



(63) 



and where and z/ ,(n ) are the coefficients of {a s /n) n in the expansion of the functions Ai 
and v l in Eq. (p2|). The two-loop expansion of the path-ordered exponential in Eq. (|50D reads 



P exp 



+ 



J J 



7T 



iln 2 (iV) 



p«i(i) x p^'C 1 ) 

5* S 



71 
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b2 ln 2 (iV)ri}5 



+ . . . 



(64) 



where is the coefficient of a s /7r in the one-loop soft anomalous dimensions of Eqs. (H3j) 

and (48). The P ordered exponential is identical to 0(a 2 s ) up to replacing T$ by T s . The 
matrices u(t 1 ,Ui,m 2 ,m 2 ) and Sij jj(1,(k, Pi, Pj) expand to one loop as 



Hij u 
Sij i j 



7T 



o(°) , " s 5(1) 

7T 



(65) 



Let us for the moment choose PIM kinematics by putting AT t = iV u = iV. We substitute 
Eqs. (|BTD, (|64"D , and ( |6~5|) into Eq. (|50D , finding the moment space expression for the NNLL- 
NNLO expanded cross section: 



(66) 



+ 



7T 



I — 

7T 



(C*' (1) + tfj (1) ) ln 2 (iV) + (Cj' (1) + Cf 1} ) ln(iV)) Tr {lliysW} 

+ Tr {4°>5W + # + (C*> (1) + ^' (1) )Tr {tff^} 

' (CI' (2) + Cf 2) + Cj (1) ^' (1) ) ln 4 (iV) 



+ (C^ + Cf (2) + 1} Cf 1} + Cj (1) Cf j ln 3 (iV) Tr {H^S^} 
- {C^ + Qp) ln 3 (iV) Tr {iff (r|'«) + sg» + //f^lf 
+ ln 2 (iV) ( C ^ + C*< a > + £ C^C**?) Tr 



fc=0 
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t „, 



-(b 2 + C\ {1) + Cf W ) ln 2 (iV) Tr (ijp) ' §® + H® 

+ (C^ + C*«) ln 2 (iV) IV {i/f^ + //^^} 
+ ln 2 (iV)Tr{i4 ) ) (I** 1 ))' (rf ^ - //^ (r™)'^ 



(o)e(o) r ii(i) r |(i) 



t~ *i y s 



With 1PI kinematics we use Eq. (p8|) and expand ln(iV u ) = In (AT) + ln(— ui/m 2 ) etc. To 
transform from moment to momentum space we use results given in the appendices of Refs. 0, 
13|, |3(J. Note that we do not keep subleading terms in our expansions when we invert to 
momentum space [EH . 



4.2 Matching 



We can identify the functions we must determine from Eq. (p6|). 

The coefficients C]^ are given in Eq. fl62|). The matrices T^j 
(|49D . The remaining unknowns are 



are given in Eqs. 



and 



S 



(0) 



ij I J i ijU 



H, 



(0) 



Tr, / P7"(°) 5(1) i rj-(l) 5(0)1 

, II T fly «Jy J 



(67) 



The leading-order soft function S'jj is simply 



S 



(o) 



tr 



(4 Cj) 



(68) 



where {c/} is our basis in color-tensor space. Recall that we have chosen the s-channel singlet- 
octet basis, Eq. (f43|). Note that = S^. The lowest-order hard function H^\j = hfXh*fj 



is calculated by projecting the lowest order amplitude Aij onto this color basis: 



h {0) 



(c K Aij^j 



l*(0) 

n ijJ 



tr [A^ c K ) (S, 



(o)V 
ii ) 



KJ 



(69) 



,(o) 



where the trace acts in ordinary color space. Note that at this stage the function h^'j may 
still depend on all other quantum numbers of the process under consideration, such as Lorentz- 
spinor and spin degrees of freedom. We are however not interested here in spin-dependent 
observables and we trace over all remaining degrees of freedom. The matrix H^\j is real and 
symmetric but may have zero eigenvalues if the Born process does not span the full color-tensor 
space, as is the case for the qq channel. For the qq channel [|TTJ the results are 



5(o) 

D qqIJ 



( N 2 

V 



o 



N 2 
— 



H, 



(o) 



qqlj 





H {0) 



(70) 



with given by 



H q%(m,m) 



2vra 2 (m) / tj + u\ 2m 2 



N 2 

c 



(71) 
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The lowest-order results are kinematics-independent. For the gg channel we find 

1) \ 

,N 2 - 4 



S 



(o) 

99 1 J 



( N C (N 2 





V 

/ rr(0) 
/ ^ggll 



H 



(0) 

gg i J 



Nr.H. 



(0) 



gg ii 



N cHg gn 
A7-2 7t(0) 

iV c -"gg 11 



2iv; 












-iV c (iV c 2 -l) J 



(72) 



(0) 



99 13 



(0) 

gg 13 



(0) 



99 13 



N r H, 



(o) 



99 13 



o-(O) , 
-"99 33 / 



(73) 



where 



Hg°g\i{m,m) 



Hg°g\ 3 {m,m) 



■na 2 {m) 



2N 2 (N- 



2^99 ' 



■na 2 s {m) 



B„ 



ti-ui 



2N C (N 2 - l) 2 99 s 2 



na 2 s {m) 



2(N, 



2-^gg 



and 



99 



It now remains to determine^] 



u\ t\ Asm 

ti U\ t\U\ 



sm 



(74) 
(75) 
(76) 

(77) 
(78) 



We do this by matching. The necessary matching conditions can be inferred by comparing the 
expansion in Eq. (|^) to exact results for the partonic cross section. 
At lowest order, the matching condition is 



Wy°(JV, s, £1, Wi, m 2 , fi 2 , a s (/j)) 



tt(0) q(0) 
n ijIJ D ijJI ■ 



(79) 



It is straightforward to check that this condition is fulfilled by our results for S^\j and H^\j. 
Neglecting 1/N terms, and choosing the case K = PIM for the moment, we have at NLO the 
matching condition 



a, ~ 



71 



f^ K ^{Hlf~S\f+H[f~Sf} = ujfj LO (N , s, ti, Mi, m 2 , /z 2 , a s (fi)) 
((C* >(1) + C|> (1) ) ln 2 (iV) + (Cj' (1) + C{ (1) ) ln(iV)) Tr {tff^} 



(80) 



71 



In(iV) Tr (4°) (rg^ + 4^«} + (Cj (1 > + «)Tr 



7 Note that we do not need the individual matrices Ssy T , and H, 
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Let us discuss this equation. To begin, the first term on the right hand side may be determined 



from the exact results in Refs. |32], |33|, [34]] , converted to moment space. To NNLL accuracy this 
term consists of all terms in the differential one-loop partonic cross sections that either diverge 
logarithmically or are constant. 

The constant terms consist, first of all, of the virtual graph contributions. In addition there 
is the soft-gluon contribution from the radiative graphs. The virtual contribution is kinematics 
independent but the soft gluon one is not because what is defined as soft differs for 1PI and 
PIM kinematics. The respective criteria are S4 < A and 1 — z < 5, with A and 5 defined in 



Eqs. and (|2"3"|). The virtual graphs contain various divergences. The ultraviolet ones are 



cancelled by renormalization of the QCD coupling constant and the heavy quark mass, for both 
kinematics in the scheme of Ref. ||36|| . Infrared divergences are cancelled by adding the virtual 
and soft contributions, so that only collinear divergences remain. These must be subtracted 
via mass factorization. The corresponding subtraction terms in general involve convolutions of 
Altarelli-Parisi splitting functions |E| with the lowest order cross section. However, we must 
mass factorize only with the soft plus virtual parts of these splitting functions. These depend on 
kinematics, see e.g. Eqs. (6.8) and (6.13) in [pUfl . What remains at NLO after mass factorization 
can now be categorized as terms multiplying ln*(A/m 2 ) or ln l (5), i — 1,2, and other terms. 
The logarithms ln(A/m 2 ) or ln(5) should be completed to plus-distributions via Eqs. (^) and 
(|23|). The final result of these procedures, in moment space, is what constitutes the first term 
on the right hand side of the equals sign in Eq. (|8~0"D . The other terms on the right hand side 
merely subtract all terms containing the singular functions in Eq. ( p2|) or (|23|). 

For 1PI kinematics there is a slight subtlety, related to the use of N t and N u in Eq. (|28"D . 
Expanding ln(A^) = ln(iV) + ln(— ti/m 2 ) (and ln(N u ) likewise) leads to extra terms containing 
ln(— ti/m 2 ) and ln(— Ui/m?) but not ln(iV) and one must be careful in accounting for such 
terms. 

We can now obtain specific results for the matching terms in momentum space. At NLO, 
T^' K on the left hand side of Eq. (^) for 1PI kinematics, 

t //ipi(Mi, «i) = ^ ( /Jpi(Mi, «i) S(s 4 ) , (81) 

„2j2_(1).S+V 



is the expression for s d a\j (s,ti,ui)/dtidui given by Eq. (4.7) of Ref. |53| for ij = qq, 



and Eq. (6.19) of Ref. |32] for ij = gg, modified by removing a factor a s /n as well as all terms 
containing ln 2 (A/m 2 ) or ln(A/m 2 ) and terms containing ln(yu/m) or ln(/i^/m)0. The definition 
of the hatted symbol in Eq. (|81|) makes the delta function explicit and is useful when we present 
the two-loop results. 

In PIM kinematics the NLO matching terms may be determined as follows. We start 
from the above 1PI results and modify those parts that are kinematics dependent. From our 
discussion following Eq. (|30"D, it is clear that these are only the soft contributions and the mass 
factorization subtraction terms. We shall label these terms by the superscript S + MF. The soft 
contributions to the PIM NLO cross sections can be computed using the results of Ref. [33. As 



mentioned in the above, we mass-factorize these results with only the soft plus virtual part of 
the splitting functions in PIM kinematics. The same procedure is done for the 1PI case using 



8 Note that in these references the factorization scale is denoted by Q and that the definition of t\ and u\ in 
Ref. |35|] for ij = qq is interchanged as compared to this work. 
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the results of Refs. [[32], RI|. For channel ij we then compute the MS quantity 



A^(M 2 ,cos#) 



2 fl),S+MF 



s, M 2 , cos6 



dM 2 dcos6 



^ S 2 
2 



dti du\ 



(82) 



PIM 



where the subscript PIM on the second term on the right hand side indicates that for t\,U\ 
one must use the expressions in Eq. ([14]) . Note that we may replace (3m by (3 to the accuracy 
at which we are working in this paper, and we shall do so for all following results. Explicit 
expressions for the terms on the right hand side of Eq. ([8^) can be found in appendix A. Then 
the PIM equivalent of Eq. (^T|) is, accounting for an overall jacobian factor originating from 
the definition of the cross section, 



7^ IM (M 2 ,cos#) 



T^ IM (M 2 , 0080)0(1-*) 



2^ipiMi>«i) 



+ Ag ) (M 2 ,cos6 



(83) 



PIM 



Note that terms containing renormalization and factorization scale dependence can be easily 
included by expanding the corresponding scale dependent exponentials in Eq. ( |50|) and adjust- 
ing the coefficients in Eq. Q6E| ) correspondingly. Alternatively, the scale dependence can be 
constructed using renormalization group methods, which we do in the next section. 

Changing to the DIS scheme at NNLL for the qq channel involves using Eq. (|53|) rather than 
Eq. (|3T| ) in the expansion, Eq. (|6T|), and extra constant terms at NLO as given in appendix B. 

Because of their length, we have collected all our results in appendix B. In the next section 
we perform a numerical study of the results obtained in this section for the inclusive partonic 
cross sections. 



5 Scaling functions 

It is convenient to express the inclusive partonic cross sections in terms of dimensionless scaling 
functions that depend only on 77, defined in Eq. (f2"U|), as 



2 



™V) = (W/i)) fe E fij' l) (v) iWAH • (84) 

m t^o £0 \ m J 

We derive LL, NLL, and NNLL approximations to fij' l \r}) for both the qq and gg channels 
using the results of the previous section, via Eqs. flD and (|15|) . Where possible we will compare 
with exact results. We also examine the results for the qq channel in the DIS scheme. Of 
course, a complete comparison between the DIS and MS schemes also requires the use of DIS 
parton densities, as discussed in the next section. 



r(fc,o) 

matics. In Fig. ^| we present the qq Born and MS NLO results for the scaling functions / ( 



We begin by showing the scaling functions (rj), k = 0, 1, 2, in both 1PI and PIM kine- 



r(0,°) 
' gg 

and fq^'" 1 , respectively. A comparison shows that, in the case of f^'"', the LL approximations 
already reproduce the exact curve quite well in both kinematics |P4| , but the NLL approxi- 



mations agree much better with the exact results to larger 77. Adding the NNLL corrections, 
numerically dominated by the negative contributions from gluon exchange between final state 
heavy quarks (Coulomb terms), leads to excellent agreement with the exact curve over a large 



range of 77. Such progressive improvement was also observed, for 1PI kinematics, in Ref. [46 
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where the absolute threshold limit was taken. We also see that, while the LL approximations 
in 1PI and PIM kinematics differ and under- or overestimate the true result, the NLL approx- 
imations in both kinematics agree better with each other, and even more so at NNLL, at least 
for 77 ^ 0.1. To summarize the one loop results in this channel, differences due to kinematics 
choice decrease as the logarithmic accuracy increases. 

The same conclusions hold for the NLO scaling function in the gg channel, fg l g '°\ shown in 
Fig. [| The Coulomb terms are positive here. At large 77, the agreement between the 1PI and 
PIM results is not as good as in the qq channel: the 1PI results at NNLL remain positive while 
the PIM ones are negative. 
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Figure 2: (a) The ^-dependence of the scaling functions ff£ (v)i fe = 0, 1 in the MS-scheme and 



1PI kinematics. We show the exact results for fj^'°\ k = 0,1 (solid lines), the LL approximation 



to fj.^ (dotted line), the NLL approximation to f^' 1 (dashed line) and the NNLL approximation 



(1.0) 



'qq 
r(l,0) 



?9 



to P q f } (dashed-dotted line), (b) The same as (a) in PIM kinematics. The spaced-dotted curve 



corresponds to the approximation involving the leading two powers of ln(/3). 



In Figs. [| and || we show the NNLO scaling functions /i| and fg 2 g '^ for both kinemat- 
ics. Again, the LL approximations differ but there is good agreement between 1PI and PIM 
kinematics at NLL accuracy. We also display the NNLL approximations for f^'°^ and f g 2 g °\ 
presently the best estimates for these functions, at least for 77 ^1, since exact calculations are 
not yet available. Near threshold, the NNLL approximations show some deviation from the 
NLL approximation due to interplay between the one-loop LL contributions and the one-loop 
Coulomb terms. Note that while the qq channel shows relatively good agreement between the 
two kinematics at large 77, the gg results are quite different since the 1PI results are positive 
while the PIM results are negative. 

We have also calculated the behavior of the scaling functions f-j'°\r]),k = 1,2, in the 
threshold limit s — > Am 2 , keeping only terms that grow as ln(/3) to next-to-leading accuracy. 
The resulting expressions are given in appendix |B.5| , and shown in Figs. 0(b), 0(b), f|(b) and 
H(b). They agree fairly well with the NLL PIM results. 

In general, comparison of the exact and approximate results in the various channels and 
kinematics in Figs. |2| and |3| indicates the range of 77 values over which we expect the threshold 
approximation to be valid. The large-77 behaviour of the scaling functions, furthest away from 
threshold, may be approximated by the methods of Ref. [47, 48, 49, 50|. However, the behavior 
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Figure 3: (a) The 77-dependence of the scaling functions fgg'°\i]), k = 0, 1 in 1PI kinematics. We 
show the exact results for fgg'°\ k = 0, 1 (solid lines), the LL approximation to fgg (dotted line), 
the NLL approximation to fgg'°^ (dashed line) and the the NNLL approximation to fgg' ^ (dashed- 
dotted line), (b) The same as (a) in PIM kinematics. The spaced-dotted curve corresponds to the 
approximation involving the leading two powers of ln(/3). 



of the scaling functions in the intermediate region, 1 ^ 77 ^ 10, can only be determined by 
exact computation. 

We now check if the qq MS results are mirrored in the DIS scheme. Since the scaling 
functions in the two schemes differ already at LL, the higher-order scaling functions will be 
quite different numerically. Of course, this difference should be compensated in principle by 
the parton densities for the scheme-independent hadronic cross section. 

Figures | and [7] show the functions and f^' ^ in the DIS scheme. We see again, as 

in Figs. [| and |4], that the exact /i* results are best approximated at small rj by the NNLL 
calculations in both kinematics while the LL and NLL calculations either over- or underestimate 
the exact calculation. The MS results generally trace the exact curves somewhat better than 
those in the DIS scheme. In particular, for the DIS scaling function fQ' , only the NNLL 
approximations provide good agreement with the exact curve for both 1PI and PIM kinematics 
over a large range in rj. This can be traced to the large delta-function terms in the scheme- 
changing functions in Eqs. ( |B.6| ) and ( |B.20| ). 

We next discuss the functions controlling the scale dependence. They can be determined 
exactly at NLO and NNLO using renormalization group methods. The exact fij'® are required 

(2 1 1 (22) 

to construct f-j ' and f^ ' . We neglect flavor mixing terms which are of order 1/N, where N is 
the moment variable. We checked that at Tevatron energies the error due to this approximation 
is less than 1% at NNLO. In this way we obtain 
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Figure 4: (a) The ^-dependence of the scaling function f^' \rf) in the MS-scheme and 1PI 
kinematics. We show the LL approximation (dotted line), the NLL approximation (dashed line) 
and the NNLL approximation (dashed-dotted line), (b) The same as (a) in PIM kinematics. The 
spaced-dotted curve corresponds to the approximation involving the leading two powers of ln(/3). 
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Figure 5: (a) The ^-dependence of the scaling function fgg'°\rf) in 1PI kinematics. We show the 
LL approximation (dotted line), the NLL approximation (dashed line) and the NNLL approximation 
(dashed-dotted line), (b) The same as (a) in PIM kinematics. The spaced-dotted curve corresponds 
to the approximation involving the leading two powers of ln(/3). 
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Figure 6: (a) The ^-dependence of the scaling function f^l'°\r]) in the DIS-scheme and 1PI 
kinematics. We show the exact result (solid line), the LL approximation (dotted line), the NLL 
approximation (dashed line) and the NNLL approximation (dashed-dotted line), (b) The same as (a) 
in PIM kinematics. 
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kinematics. We show the LL approximation (dotted line), the NLL approximation (dashed line) and 
the NNLL approximation (dashed-dotted line), (b) The same as (a) in PIM kinematics. 
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where Pgg\ Pjj$ and P}$, are the one- and two-loop splitting functions [41], |52| , and 
~ indicates the neglect of flavor-mixing terms. The convolutions involving a scaling function 
fi*' ^ are defined as 



(/if ® P$) (TJ{X)) = [ dz ( V (XZ)) P<p k (z) 



(91) 



4m 2 /s 



with rj(x) = xs/(Am 2 ) — 1. The standard convolution of two splitting functions, P^' ' and 
in Eqs. © and © are 



,(0,0), 



X = 



with i 
Ref, 



naturally agree with the results for 



(92) 



and fol'^ in 



Equations (H) and vt ^, _ ^ ...... ,,,,, 

31]. We have also checked that the above results agree to NNLL with the expressions 
in appendix B when integrated as in Eqs. (§) and flTsp. Note that in appendix B we have 
also given results for terms in the differential NNLO functions controlling the scale dependence 
beyond NNLL accuracy, thus deriving all the soft plus virtual terms in these functions. The 
exact hard terms are calculated only for the integrated cross section as above. 

We begin by showing the scale-changing scaling functions in the qq channel and MS scheme, 
comparing 1PI and PIM kinematics. In Fig. || we show and note that the PIM LL 

approximation reproduces the exact curve somewhat better than the 1PI LL aproximation. 
The NLL approximations agree better, even for larger rj. 

In Fig. H we show the NNLO scaling functions f^' 1 ^ and ■ We compare the exact curves 
calculated from Eqs. 



and (37) with our LL, NLL, and NNLL approximations. Again we 
see that the NNLL approximations provide a remarkably good description of the exact results, 
both in shape and magnitude. The NNLL curves for 1PI and PIM kinematics are in very good 
agreement with each other, i.e. ambiguities from the kinematics choice are very mild. Similar 
conclusions hold for the gg scaling functions fjjl' 1 ', fgg' 1 ^, and fg 2 g ' 2 \ shown in Figs. [10] and [11. 
Note however that also here the NNLL 1PI and PIM results for fgg' 1 ^ differ at large r\. The 
PIM NNLL scaling function differs significantly from the exact result. 

( 2 1) f 2 2) 

For completeness, we also display the NNLO scaling functions /V and f^' in the DIS 
scheme. The NLL approximations roughly trace the exact curve. The (percentage-wise) large 
difference between the NLL approximation and the exact curve close to threshold may be 
attributed to the large constants in the one-loop scheme-changing functions in Eqs. ( P-6| ) and 
(|B.20j ) that interfere with the one-loop LL terms. These are accounted for at NNLL accuracy, 

(2 1) 

as Fig. [T2| demonstrates. However, the NNLL results for f^' differ more significantly between 
the two kinematic choices than their MS counterparts. 
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Figure 8: (a) The ^-dependence of the scaling function f^'^^rj) in 1PI kinematics. We show the 
exact result (solid line) , the LL approximation (dotted line) and the NLL approximation (dashed line) . 
(b) The same as (a) in PIM kinematics. 
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Figure 9: (a) The r/-dependence of the scaling functions f^' (ry) (MS-scheme) and f^' (rj) in 1PI 
kinematics. We show the exact results (solid lines), the LL approximations (dotted lines), the NLL 
approximations (dashed lines) and the NNLL approximations (dashed-dotted lines), (b) The same as 
(a) in PIM kinematics. 
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Figure 10: (a) The jy-dependence of the scaling function fgl'^irj) in 1PI kinematics. We show 
the exact result (solid line), the LL approximation (dotted line) and the NLL approximation (dashed 
line), (b) The same as (a) in PIM kinematics; LL and NLL approximation coincide here. 
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Figure 11: (a) The ^-dependence of the scaling function fg g ' (rj), I = 1,2 in 1PI kinematics. We 
show the exact results (solid lines), the LL approximations (dotted lines), the NLL approximations 
(dashed lines) and the NNLL approximations (dashed-dotted lines), (b) The same as (a) in PIM 
kinematics. 
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Figure 12: (a) The //-dependence of the scaling function f\J '(r/) in the DIS-scheme and IPI 
kinematics. We show the exact result (solid line), the LL approximation (dotted line), the NLL 
approximation (dashed line) and the NNLL approximation (dashed-dotted line), (b) The same as (a) 
in PIM kinematics. 



6 Hadronic cross sections 

In the previous section we examined the quality of our threshold approximations at the parton 
level. Here we assess these approximations at the hadron level for inclusive top quark pro- 
duction at the Fermilab Tevatron and bottom quark production at HERA-B. The inclusive 
hadronic cross section is the convolution, Eq. (|i~9]) , of parton distribution functions with the 
partonic cross section, expressed in terms of the scaling functions, Eq. (Q). To facilitate the 
understanding of the results in this section in terms of those of the previous section, we plot the 
flux factors $^(77, // 2 ), Eq. (]T8|), for the above cases in Fig. [13]. They show which rj values re- 
ceive the most weight in the convolution integral. In our numerical studies we use the two-loop 
expression of a s and the CTEQ5M (MS scheme) or CTEQ5D (DIS scheme) parametrizations 
of the parton distributions |k| not only for the NLO results, but also for the NNLO (NNLO 
parton distributions are not yet available) and LO results. Thus in this section we keep the 
nonperturbative part of our results fixed when studying the effect of increasing the perturbative 
order of our partonic cross sections. For top quark production at the Tevatron and bottom 
quark production at HERA-B the calculations probe the moderate to large x region where the 
parton distributions are well known, see Fig. [13]. We use five and four active flavors respec- 
tively for these cases and fix fiR = \i. Except where specified otherwise we have multiplied 
the non-exact scaling functions at NLO and NNLO with a damping factor 1/^/1 + rj, as in 
Ref. f54| in order to lessen the influence of the large i] region of the scaling functions where 
threshold logarithms become less dominant and we lose some theoretical controf]. Figure [T3 



demonstrates that this factor indeed damps the large rj region of the NNLO scaling functions 
fij'°\ while leaving the small and medium r/ regions unaffected. The effect of the damping 
factor will also be made explicit in the tables. 



9 In IPI kinematics such a factor is effectively equivalent to including 9(s4 — m 2 ) in Eq. (g), see Ref. 
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Figure 13: (a) The qq parton flux factor ln(10)ry/(l + r?)$^ (see Eq. (gJD) for the CTEQ5M 
parametrization at the Tevatron (upper three curves, y/S = 1.8 TeV and m = 175 GeV) and for 
HERA-B (lower three curves \/~S = 41.6 GeV and m = 4.75 GeV). We show results for fi = m (solid 
curves), = m/2 (dotted curves), and fi = 2m (dashed curves), (b) Same as (a) for the gg parton 
flux factor ln(10)r// (1 + r})$> gg . Now the upper set of curves correspond to HERA-B and the lower set 
to the Tevatron. 
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without the damping factor 1/^/1 + rj. We show the NNLL approximation to f^'°^ in 1PI/PIM kine- 
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6.1 Results for tt production at the Tevatron 

We first discuss top quark production at the Tevatron in proton-antiproton collisions. (For a 
review of the data see Refs. ]5R We give results for both Run I (y/S = 1.8 TeV) and Run 
II (y/S = 2.0 TeV). At the Tevatron, top quarks are mainly produced in pairs. With a top 
quark mass of 175 GeV, the dominant production channel is qq annihilation, constituting about 
90% of the total cross section at the Born level for y/S =1.8 TeV, with the gg channel making 
up the remainder. The NLO corrections in the qq channel are moderate, of the order of 20%, 
whereas those in the gg channel are more than 80% so that the gg channel gains significance at 
NLO0. These large corrections originate predominantly from the threshold region. Figure |L3 
shows that the range 0.1 ^ r\ ^ 2 contributes the most. At still higher orders, this trend 
continues: the relative corrections to the gg channel are larger than those for the qq channel, 
as we shall see. This is due to the larger color factors in the analytical expressions for the 
corrections for the gg channel. As we mentioned earlier, the qg and qg channels give negligible 
contributions to the total cross section and are not considered here. 

We begin by comparing MS results in 1PI and PIM kinematics, including only the qq channel 
in Eq. fll9|), at v^S* = 1.8 TeV. In Fig. [L5] we show the Born cross section and the exact and 



approximate NLO corrections, the latter at both NLL and NNLL accuracy, for fx — m in the 
range 150 < m < 200 GeV. We see that our NLO 1PI approximations are a little larger than 
the exact result while the NLO-NNLL approximation in PIM kinematics is indistinguishable 



from the exact answer. The corrections are about 20-30% in the mass range shown. In Fig. 16 
we give the equivalent results for the gg channel. Here the 1PI approximations agree with the 
exact result better than the PIM ones. 

Figure |17|(a) displays the approximate qq NNLO corrections at NLL and NNLL accuracy for 
jj, = m as a function of the top quark mass in a direct comparison of the 1PI and PIM results. 
The shaded area indicates the kinematics ambiguity at NNLL, about 0.6 pb at m = 175 GeV. 
The figure shows that the NNLL ambiguity is larger than the NLL one. (Thus the small size of 
the NLL kinematics ambiguity seems somewhat accidental.) Figure |l7|(b) shows similar results 
for the gg channel. 

For completeness, we show the corresponding results for the qq channel in the DIS scheme in 
Figs. [T3| and |TU|. Again PIM kinematics approximates the exact results somewhat better than 
1PI kinematics at NLO-NNLL. We see that the NNLL-NNLO kinematics ambiguity in Fig. 0, 
again indicated by the shaded region, is greatly reduced compared to the MS case, Fig. |T7|(a). 
It should however be kept in mind that, in the DIS scheme, the parton densities absorb large 
threshold logarithms which are not properly accounted for at NNLO if one uses NLO parton 



distributions as in Fig. [19]. Therefore it seems likely that the DIS scheme kinematics ambiguity 
is somewhat underestimated in Fig. [19|. 

In Fig. ^ the sum of the qq and gg channels in the MS scheme is shown as a function of 
the top quark mass. We display the exact NLO cross section and the approximate NNLO cross 
section, which is the sum of the exact NLO cross section and the NNLL-NNLO corrections. We 
show results for \x = m/2,m, and 2m. We see that the NNLO cross section is uniformly larger 
than the exact NLO one, although less so in PIM kinematics, and that the scale dependence of 
the NNLO cross section is considerably reduced relative to NLO. Comparing Figs, ^(a) and 
(b) we observe that the kinematics ambiguity is larger than the scale uncertainty. 

We now turn to a more detailed study of the scale dependence of the inclusive MS top 
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Recall that we also use NLO parton distribution function for LO cross sections. 
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Figure 15: (a) The qq channel contribution to the top quark cross section at the Tevatron with 
y/S = 1.8 TeV and ji = m in the MS scheme. We show the Born term (upper solid line), the exact NLO 
corrections (lower solid line) and the 1PI approximate NLL (dashed line) and NNLL (dashed-dotted 
line) one-loop corrections, (b) The same as (a) in PIM kinematics. 
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Figure 16: (a) The gg channel contribution to the top quark cross section at the Tevatron with 
= 1.8 TeV and \x = mm the MS scheme. We show the Born term (upper solid line), the exact NLO 
corrections (lower solid line) and the 1PI approximate NLL (dashed line) and NNLL (dashed-dotted 
line) one-loop corrections, (b) The same as (a) in PIM kinematics. 
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Figure 17: (a) The qq channel two-loop corrections to the top quark cross section at the Tevatron 
with y/S = l.S TeV and H = m in the MS scheme. We show the 1PI approximate NLL (upper dashed 
line) and NNLL (upper boundary of shaded region) two-loop corrections and the PIM approximate 
NLL (lower dashed line) and NNLL (lower boundary of shaded region) two-loop corrections, (b) The 
same as (a) for the gg channel two-loop corrections. 
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Figure 18: (a) The qq channel contribution to the top quark cross section at the Tevatron with 
= 1.8 TeV and \x = m in the DIS scheme. We show the Born term (upper solid line), the exact NLO 
corrections (lower solid line) and the 1PI approximate NLL (dashed line) and NNLL (dashed-dotted 
line) one-loop corrections, (b) The same as (a) in PIM kinematics. 
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Figure 19: The qq channel contribution to the top quark cross section at the Tevatron with 
VS = 1.8 TeV and [i = m in the DIS scheme. We show the 1PI approximate NLL (upper dashed line) 
and NNLL (lower boundary of shaded region) two-loop corrections and the PIM approximate NLL 
(lower dashed line) and NNLL (upper boundary of shaded region) two-loop corrections. 
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Figure 20: (a) The top quark cross section at the Tevatron with v 7 *? =1.8 TeV for the sum of the 
qq and gg channels in the MS scheme. We show the exact NLO cross section for fi = m (solid line), 
m/2 (upper dotted line), and 2m (lower dotted line), and the 1PI approximate NNLL-NNLO cross 
section for \x = m (dashed-dotted line), m/2 (upper dashed line), and 2m (lower dashed line), (b) 
The same as (a) in PIM kinematics. 
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cross section. We show the sum of the qq and gg channels in Fig. 21 at v^S" =1.8 TeV and 
m =175 GeV at several orders and to different accuracies. The Born and NLO results shown 
are exact. One of the NNLO curves is constructed by adding the contributions from the NNLL 

(2,0 



approximate two- loop scaling functions , I 

by adding instead the exact f^' 1 ^ and fj; 2 ' 2 ^ functions to the approximate f-j'°\ 



0,1,2, to the exact NLO results, the other 

thus making 

the changes to the partonic cross section when changing [i exact. The differences between 
these two NNLO curves are due to subleading terms and represent, for each kinematics, the 
corresponding ambiguity in the scale dependence. At very small /i the contributions of the 
terms involving /£• and fij' 2 ^ are much larger that the contribution from fj? . The sizable 
difference between the two NNLO curves in Fig. ^l](a) is in fact mainly due to the difference 
between the exact and NNLL 1PI results for jfg| at medium and large r] in Fig. |9|(a). Even if 

(2 1) 

we include all soft plus virtual terms in the approximate 1PI , as derived in appendix B, 
there is still a sizeable difference from the exact result. Therefore this difference stems mainly 

(2 1) 

from the hard, i.e. 0(1/ N), terms in / g y . 

The NNLO differences in PIM kinematics are much smaller than in 1PI kinematics, in 

(2 1) 

correspondence with the good agreement of the exact and NNLL PIM results for /V over all 



relevant rj, as shown in Fig. f21|(b). 
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Figure 21: (a) The /i-dependence of the top quark cross section at the Tevatron with y/S = 1.8 
TeV and m = 175 GeV for the sum of the qq and gg channels in the MS scheme. We show the Born 
(upper solid line at small n/m) and the exact NLO (lower solid line at small fi/m) cross sections, the 
1PI approximate NNLL-NNLO cross section (dashed line) and the NNLO estimate with only ffQ 
and fgg'°^ NNLL approximate (dashed-dotted line), (b) The same as (a) in PIM kinematics. 



Next, in view of the upgrade in energy for the Tevatron from \/~S = 1.8 to 2.0 TeV, we 



investigate the v S dependence of the top quark production cross section. In Fig. |22] we present 
the inclusive cross section for the sum of the qq and gg channels in the MS scheme at NLO 
and NNLO as a function of a/5*. The NNLO curves result from adding the NNLL NNLO 
corrections to the exact NLO cross section. We have normalized all calculations to the value 
of the exact NLO cross section at fi = m. Comparing 1PI and PIM kinematics we find that at 
lower energies, where the qq channel is dominant, the NNLO results are 10-30% larger than at 
NLO in both kinematics. As \/S increases, the gg channel, with its larger corrections, grows 
in importance. The cross sections are also more sensitive to the large r] behavior of the scaling 
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functions. This leads to large kinematics differences at energies above 5 TeV although the scale 
uncertainties remain small. In PIM kinematics this is due to fig' ' which is large and negative 
for 7] > 0.1 at NNLL, thus reducing the NNLO results relative to the exact NLO cross section. 
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Figure 22: (a) The top quark cross section with m = 175 GeV at pp colliders for the sum of the qq 
and gg channels in the MS scheme as a function of normalized to the exact NLO cross section at 
fj, = m. We show the ratios of the exact NLO cross sections for [i = m (solid line), m/2 (upper dotted 
line), and 2m (lower dotted line), and the ratios of the 1PI approximate NNLL-NNLO cross sections 
for /i = m (dashed-dotted line), m/2 (upper dashed line), and 2m (lower dashed line), (b) The same 
as (a) in PIM kinematics. 



We now present the values of the NLO and NNLO total MS cross sections for top quark 
production at the Tevatron with v^S* = 1.8 TeV and 2.0 TeV for m = 175 GeV and \i = m, m/2, 
and 2m in Tables [I] and 0. They detail the effects of both multiplying the approximate scaling 
functions fjj , k = 0, 1,2, with the damping factor 1 / y/1 + rj and using the approximate or 

exact (Eqs. flgop-flSOD) scaling functions /i 2 and f^' 2 ^ on the NNLO results. Note that all 
NNLO results in the figures correspond to the option AP+DF, except in Fig. |21] where we also 
show NNLO results corresponding to the EX+DF choice (dash-dotted lines). Comparing the 
NNLO and NLO cross sections in the tables shows that the NNLO corrections are much smaller 
in PIM kinematics than in 1PI kinematics, as already shown in Figs, ^(a) and (b). We see 
that for the top cross section the effect of the damping factor is rather small. We also note that 
the use of exact f^' 1 ^ and fjj' 2 ^ affects the \i = m/2 case the most. The difference between 
the exact and approximate calculation is larger for 1PI than PIM kinematics, as can also be 



observed in Fig. 21 



We note that AP results in 1PI kinematics were also presented in Ref . [[HJ . There are some 
small numerical differences with the results in this paper stemming mainly from using slightly 
different analytical expressions, all equivalent at threshold. To be specific, the expressions 
t\ + u\ and s 2 — 2tiU\ are equivalent at threshold and either choice can be made in our NNLO 
expansions. Different choices produce slightly different numerical results away from threshold 
and the variation in these results represents a small but inherent uncertainty in the cross 
sections. 

Let us finally comment on the applicability of our results for top quark production at the 
LHC where = 14 TeV. At this pp collider the gg channel is dominant (about 90% of the 
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y/S = l. 


S TeV 






\i = m/2 


/i = m 


ji = 2m 


NLO 


5.39 


5.16 


4.61 


NNLO 


;ipi) 


AP 


6.27 


6.32 


6.12 


AP+DF 


6.19 


6.22 


5.97 


EX 


6.91 


6.32 


5.84 


EX+DF 


6.73 


6.22 


5.78 


NNLO (PIM) 


AP 


5.37 


5.35 


5.26 


AP+DF 


5.45 


5.45 


5.29 


EX 


5.18 


5.35 


5.26 


EX+DF 


5.36 


5.45 


5.31 



Table 1: The hadronic it production cross sections in pb for pp collisions in the MS scheme with 
= 1.8 TeV and m = 175 GeV, for fx = m,m/2, and 2m. The labelling of the NNLO results 
corresponds to with or without the damping factor (DF) + rj, and using exact (EX) or NNLL 

approximate (AP) scaling functions f^' 1 ^ and f^ 2 ' 2 ^ . 



y/S = 2.0 TeV 




H = m/2 


\i = m 


fi = 2m 


NLO 


7.37 


7.10 


6.36 


NNLO (1PI) 


AP 


8.58 


8.71 


8.46 


AP+DF 


8.47 


8.56 


8.24 


EX 


9.47 


8.71 


8.06 


EX+DF 


9.21 


8.56 


7.97 


NNLO (PIM) 


AP 


7.27 


7.27 


7.17 


AP+DF 


7.40 


7.43 


7.24 


EX 


6.88 


7.27 


7.20 


EX+DF 


7.20 


7.43 


7.29 



Table 2: The hadronic tt production cross sections in pb for pp collisions in the MS scheme with 
y/S = 2.0 TeV and m = 175 GeV, for fx = m,m/2, and 2m. The labelling of the NNLO results 
corresponds to with or without the damping factor (DF) l/\/l + rj, and using exact (EX) or NNLL 
approximate (AP) scaling functions fj?' 1 ^ and f^' 2, '. 
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total) because only sea quarks contribute to the antiquark distributions in the qq channel. 
Although top production at these energies is far from the hadronic threshold region, it might 
be close enough to partonic threshold for threshold resummation to be relevant since the gluon 
flux may favor small values of s = X1X2S <C S (see the discussion in Refs. [^7| |58j). Fig. ^3 
indeed confirms this, but also shows that the hadronic cross section is sensitive to the high- 
energy behaviour of the scaling functions. Thus, estimates of the inclusive top cross section at 
the LHC based on the threshold approximation alone are unreliable. 
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Figure 23: (a) The qq parton flux factor, ln(10)r//(l + r])(& q q, for the CTEQ5M parametrization at 
the LHC (VS = 14 TeV and m = 175 GeV). We show results for /i = m (solid curve), fi = m/2 (dotted 
curve), and \i = 2m (dashed curve), (b) Same as (a) for the gg parton flux factor, ln(10)r//(l + rj)$> gg . 



6.2 Results for bb production at HERA-B 

In this section we present the inclusive cross section for bottom quark production at fixed- 
target pp experiments, in particular, the HERA-B experiment. The energy of the proton beam 
at HERA-B is 920 GeV so that = 41.6 GeV. Here the gg channel is dominant (about 70% 
of the total cross section), with the qq channel contributing the remainder. The gq and gq 
channels are again negligible, of the order of a few percent. Figure [13] shows that the r] ^ 1 
region is dominant in the convolution with the parton densities, Eq. (|19"D. 

In Fig. |2| we present our NLO and NNLL-NNLO results for the b quark production cross 



section in fixed-target pp interactions as a function of beam energy in the range 200-1200 
GeV/c with m = 4.75 GeV. Comparing 1PI and PIM kinematics we find that, particularly at 
high energies, the NNLO predictions are different. Again, this can be attributed to increased 
sensitivity to the high-energy asymptotics of the scaling functions. In Table ^]we list the NLO 
and NNLO inclusive b quark production cross sections. We show the effects of both using the 
damping factor 1/yT+ly and the approximate or exact (Eqs. (|85|) -(90)) scaling functions f^f' 1 ^ 

and fjj' 2 ^ on the NNLO results. 

We observe that the damping factor has a significant effect at small scales. The use of exact 
scaling functions also has a dramatic effect since it causes the NNLO bottom cross section in 
PIM kinematics to become negative at /i = m/2. The 1PI inclusive cross section at NNLO is 
significantly larger than the PIM cross section. Only the PIM NNLO results show a significant 
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Figure 24: (a) The total MS bottom quark cross section at fixed-target pp experiments with 
m = 4.75 GeV. We show the exact NLO cross section for fi = m (solid line), m/2 (upper dotted line), 
and 2m (lower dotted line), and the NNLL-NNLO cross section for fi = m (dashed-dotted line), m/2 
(upper dashed line), and 2m (lower dashed line), (b) The same as (a) in PIM kinematics. 



^ = 41.6 GeV 




fi = m/2 


fi = m 


fi = 2m 


NLO 


31.1 


17.8 


9.7 


NNLO 


;ipi) 


AP 
AP+DF 

EX 
EX+DF 


77.1 
72.2 
91.1 
80.7 


39.1 
37.2 
39.1 
37.2 


21.9 
20.8 
20.8 
20.3 


NNLO (PIM) 


AP 
AP+DF 

EX 
EX+DF 


18.1 
27.3 
-19.8 
-2.5 


19.0 
21.8 
19.0 
21.8 


14.0 
14.7 
15.6 
16.2 



Table 3: The hadronic bb production cross sections in nb at HERA-B in the MS scheme with y/S = 41.6 
GeV and m = 4.75 GeV for fi = m,m/2, and 2m. The labelling of the NNLO results corresponds 
to with or without the damping factor (DF) 1 / y/l + rj, and using exact (EX) or NNLL approximate 
(AP) scaling functions f^' 1 ^ and fff' 2 ^- 
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reduction in scale dependence. Note also that for most options the scale uncertainty is larger 
than the kinematics ambiguity. In general, we see that the theoretical control over this cross 
section is still rather poor. 



7 Conclusions 

In this paper we have resummed threshold enhancements to double-differential heavy quark 
hadroproduction cross sections, both in single-particle inclusive and pair-invariant mass kine- 
matics and in the qq and gg production channels. We used these resummed expressions to 
construct analytic approximations for the heavy quark cross sections up to NNLO, including 
the first three powers (NNLL) of the large logarithms. This involved including color-coherence 
effects and soft radiation attached to all the one-loop virtual corrections, the contribution of 
which we determined via matching conditions. We derived exact results for the scale-dependent 
terms at NNLO using renormalization group methods. 

We examined the magnitude of the NNLO correction to the inclusive cross section at the 
parton and hadron levels and studied its variations due to kinematics choice, factorization 
scheme choice, and changes in the renormalization/factorization scale for top quark production 
at the Tevatron and bottom quark production at HERA-B. We found that the scale dependence 
of the cross section after including the NNLO correction is typically substantially reduced, as 



expected on general grounds [|9|, |6( 

We now provide what are, in our judgement, reasonable values of the approximate NNLO 
top cross section at 1.8 and 2.0 TeV and the approximate NNLO bottom cross section at 41.6 
GeV. These values, based on a combination of the results and their uncertainties, are obtained 
in the following procedure. For each energy, we take the AP+DF results for each kinematics 
choice, given in the tables in section 6. The scale uncertainty for 1PI or PIM is given by the 
maximum and minimum values of the cross sections in each row. The central value of the cross 
section is the average of the /i = m values in 1PI and PIM kinematics. Two errors are assigned 
to the resulting cross section. The first is due to the kinematics-induced ambiguity and is the 
difference between the central value and the /i = m values obtained in 1PI and PIM kinematics 
alone. The second is the weighted average of the scale uncertainties in the two kinematics, 
giving more weight to smaller scale uncertainties. Thus we obtain the following NNLO top 
production cross sections at the Tevatron, 

0- tt -(1.8TeV) = 5.8 ±0.4 ±0.1 pb , (93) 

and 

<x tt -(2.0TeV) =8.0 ±0.6 ±0.1 pb . (94) 
The HERA-B cross section to NNLO is 

0-65(41.6 GeV) = 30 ± 8 ± 10 nb . (95) 

Recall that the first set of errors indicates the kinematics ambiguity while the second is an 
estimate of the scale uncertainty. Note that the scale uncertainty is considerably smaller than 
the kinematics uncertainty for top production. There are other sources of uncertainties in 
the total cross section such as parton distribution function uncertainties |62|, |63], [>|, |65 



ambiguities in the analytical expressions near threshold, as discussed in the previous section, 
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and contributions from subleading logarithms pOl. Contributions from even higher orders may 



also be significant |30|. Due to the relative large uncertainties in the cross section we give one 
less significant figure in our combined results than in the tables of the previous section. 

We note that the central values for the cross sections do not change at all for top production 
and change only slightly for bottom production if we use any of the other options in the 
tables. The kinematics uncertainty also changes only slightly for the other options. The scale 
uncertainty is increased substantially, however, if we use the EX or EX+DF options. The 
combined results presented here reflect our subjective judgment - more unbiased results are 
presented in section 6. 

In conclusion, we have provided, using soft gluon resummation techniques, improved finite- 
order perturbative estimates for heavy quark production along with their uncertainties. We 
hope that our results will thereby allow more meaningful comparisons with measurements. 
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A NNLL matching terms for PIM kinematics 

Here we collect the terms required for NLO matching in PIM kinematics. We first define the 
Born functions that occur in these terms. Recall that in our notation the superscript (n) 
indicates the power, (a s (fj,ji) / 7r) n , in the expansion of the corresponding quantity. 
The qq channel 1PI Born function is 

F^is^m) = Traft*)^ (^±^+ 2 Jpj . (A.l) 

The PIM Born function is 

F| om ( S ,M 2 ,cos#) = ^F| OTn ( S ,t 1)Ml )l PIM , (A.2) 

where the subscript PIM indicates that the expressions for t\ and U\ in Eq. flTjj) should be used 
with p M replaced by (5. 

The gg channel 1PI Born function is 



Mi ti Asm 2 ( sm 2 \ 
B qq = — + — + 1 . (A.4) 
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The color factors Cq,Ck and, for later use, Cqed, are defined as 



N 2 - 1 

Co = N C (N 2 - 1) , C K =■ ' ' 



Nr. 



N 4 - 1 

Cr 



y QED 



AT2 



The PIM Born function is 

F„ B , orn (s,M 2 ,cos6 



fC rn ( S »*l»«l)|P!M 



(A.5) 



(A.6) 



We again indicate that the expressions in Eq. (0) should be used for t\ and U\ with /3 M replaced 
by (3. We also define for PIM kinematics 



BL = £b. 



yy 



yy 



(A.7) 



Let us now turn to the one-loop S + MF cross sections, as defined in section |]. The MS 
scheme 1PI qq result is 
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with x = (1 - /3)/(l + /?). The MS scheme PIM gg result is 
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The MS scheme 1PI gg result is 
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The MS scheme PIM result is 
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B NNLL differential partonic heavy quark cross sections 
to NNLO 



In this appendix we collect all our results for the NNLL-NNLO differential partonic heavy 
quark cross sections. We first make some general observations. The NLO 1PI results agree 
with those in Refs. |33], The NLO PIM results are consistent with Ref. ||34j| . The qq 
channel results contain terms that are antisymmetric under ti,ui interchange, a consequence 
of charge conjugation asymmetry of the initial state. The gg channel results are symmetric 
under t\ <-> u\. In the qq channel the Born function factors out as a whole for LL and NLL 
terms. In the gg channel the Born function function factors out as a whole only for LL terms. 
We also use the notation that the superscript (n) denotes the coefficient of (o> s (hr) /ii) n in the 
expansion of the corresponding quantity. 



B.l The qq channel in 1PI kinematics 

In the qq channel, Eq. (HI), the Born cross section is 
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with F^- orn (s, ti, Mi) given in Eq. ( |A.1|) . The MS one-loop NNLL corrections are 
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The MS two-loop NLL corrections are 
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To achieve NNLL accuracy one must add 
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where we suppressed the channel label on the soft anomalous dimension matrix elements. 

We note that at NNLL accuracy we have derived all NNLO soft plus virtual terms involving 
the scale, except for <5(s 4 ) terms involving single logarithms of the scale. We can go beyond 
NNLL accuracy and also derive these terms in the partonic cross section by requiring that the 
scale dependence in the hadronic cross section cancel out. These terms are 
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with T ,( lg 1P1 (s,ti,Ui) = T^\ Pl (s,ti,Ui) / F^ orn (s,t 1 ,ui). We have checked that these results are 
consistent with the exact expressions in Eq. ( |86D and with the expansion of the resummed cross 
section beyond NNLL accuracy as discussed in Ref. J3D . 
The DIS one-loop NNLL corrections are 
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The DIS two-loop NLL corrections are 
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To achieve NNLL accuracy for the DIS scheme at two loops we have to add to the NNLL MS 
terms 
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The gg channel in 1PI kinematics 

In the gg channel, Eq. (|5]), the Born cross section is 
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The MS two-loop NLL corrections are 
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To achieve NNLL accuracy one must add 
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where we suppressed the gg channel label on the soft anomalous dimension matrix elements. 

As we discussed in the previous subsection, we can also derive the NNLO S(si) terms 
involving single logarithms of the scale. These terms are 
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with T' gglPl (s, t\, Mi) = T^- ) 1PI (s, ti, -ui)/F^ orn (s, ti, Mi). We have checked that these results are 
consistent with the exact expressions in Eq. fl39| ) and with the expansion of the resummed cross 
section beyond NNLL accuracy. 



B.3 The qq channel in PIM kinematics 

In the qq channel, Eq. fll2|) , the Born cross section is 



d 2 a^(s,M 2 ,cos6) p 



,(°)/ 



dM*dcos9 2^ (S '* 1)U1)|PIM 

with F| orn (s,M 2 ,cos#) defined in Eq. ( gg) . 
The MS one-loop NNLL corrections are 

- 99 V 7 -^i ) (l-^,s,M 2 ,cos6 



(iM 2 cicos# 

= -F| om ( S ,M 2 ,cos£)) 



X 



s 

4C F 



ln(l - z) 
1-z 



-5(1 - z) F* oin (s, M 2 ,cos6 
s 



P (l) / . \ | 

2^99 (S4,S,ti,tti)|piM 



(B.15) 



(B.16) 



+ 



■ 1 " 


.1 


[l-z. 





-31n(^ 



In 



sm 



+ ReL 



47 
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The MS two-loop NLL corrections are 
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To achieve NNLL accuracy to two-loops one must add 
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As for 1PI kinematics, we note that at NNLL accuracy we have derived all NNLO soft plus 
virtual terms involving the scale, except for 5(1 — z) terms involving single logarithms of the 
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scale. We can go beyond NNLL accuracy and derive these terms in the partonic cross section 
by requiring that the scale dependence in the hadronic cross section cancel out. These terms 
are 
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with T'^ PIM (M 2 ,cos#) = sf^ ) PlM (M 2 ,cos6)/F^ orn (s,M 2 ,cos6). Again, we have checked that 
these results are consistent with Eq. (|86|) and with the expansion of the resummed cross section 
beyond NNLL accuracy. 

The DIS one-loop NNLL corrections are 
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To achieve NNLL accuracy for the DIS scheme at two loops we have to add to the NNLL MS 
terms 



d 2 a ( 9{s,M 2 , cost 



NNLL terms 



d 2 a ( 9{s,M 2 , cost 



DIS 



dM 2 dcos6 

-2C F f£ Fm (M 2 ,cos9) 
"ln(l - z 



dM 2 dcosO 



NNLL terms 



MS 





"ln(l -z) 


MS 


1 - z 



+ -F* oin (s,M 2 ,cost 



1 - z 



l - C F K + 6C F (ReT^ -C F + C F In (2^/^ 



(B.22) 



49 



3 
-( 

2 



l-z 



m 
45 



- ^C F b 2 + 2C F b 2 In ( — ) + 16C 2 C 2 + -^C 2 



B.4 The channel in PIM kinematics 

In the gg channel, Eq. flilf), the Born cross section is 
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with F^ orn (s, M 2 , cos^) defined in Eq. ( |A.6| ). The MS one-loop NNLL corrections are 
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The MS two-loop NLL corrections are 
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As we discussed in the previous subsection, we can also derive the NNLO 8(1 — z) terms 
involving single logarithms of the scale. These terms are 
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with T'^' PIM (M 2 ,cos#) = sf; ( J ) p IM (M 2 ,cos^)/F Born (s,M 2 ,cos^). Again, we have checked that 
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these results are consistent with Eq. (|8^) and with the expansion of the resummed cross section 
beyond NNLL accuracy. 



B.5 Results for s —> Am 2 

Here we present the scaling functions of section |5| up to two loops for the inclusive cross section 
for s — > 4m 2 . They may be obtained from the results in sections B.l to B.4 via Eqs. (||]) and 
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(15), keeping only terms that behave as ln(/3). We checked that we find the same results for 
both kinematics. Our results are as follows. 
For the qq channel in the MS scheme 
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while in the DIS scheme we have 
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For the gg channel we find 
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